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Abstract 

Recently the BES collaboration has announced observation of a resonant state in the Tr + ir~ri' 
spectrum in J/ip — ► 77r + 7r - ?/ decay. Fitting the data with a ^ state, the mass is determined to 
be 1833.7 MeV with 7.7a statistic significance. This state is consistent with the one extracted from 
previously reported pp threshold enhancement data in J ftp — * 'ypp. We study the properties of this 
state using QCD anomaly and QCD sum rules assuming X(1835) to be a pseudoscalar and show 
that it is consistent with data. We find that this state has a sizeable matrix element < 0\GG\G p > 
leading to branching ratios of (2.61 ~ 7.37) x 10" 3 and (2.21 ~ 10.61) x 10~ 2 for J/ip -> ~/G p and 
for G p — ► ir + ir~rj', respectively. Combining the calculated branching ratio of J/ip — > ^G p and data 
on threshold enhancement in J/ip — ► jpp, we determine the coupling for G p — p — p interaction. 
We finally study branching ratios of other J/ip — > 7 + three mesons decay modes. We find that 
J/ip — ► jG p — > 7(7T + 7r _ ?7, KK-k ) can provide useful tests for the mechanism proposed. 
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Recently the BES collaboration has announced observation of a resonant state in the 
TT + -K~r]' spectrum|l| in J/ip — > / yir + -K~ri'. A fit for a 0~ + resonant state with a Breit- 
Wigner function yields a mass M = 1833.7 ± Q.5(stat) ± 2.7{syst) MeV, a width T = 
67.7 ± 20.3(stat) ± 7.7(syst) MeV and a production branching ratio B(J/ip — > 7X)5(X — > 
7T + 7T~r]') = (2.2 ± 0.4(stat) ± 0.4(s?/si)) x with a 7.7<r statistic significance. The mass 
of this state is consistent with that extracted from the enhanced threshold pp events inQ] 
J/ip —>■ jpp- The properties of this state cannot be explained by known particles. Various 
models have been proposed to explain the possible resonance P, U, ^]. Further experimental 
confirmation of this state is needed. 

There are some hypothetic candidate particles which may fit in the picture. Some of 
the possibilities include a pp bound statejj], and a pseudoscalar glueball statej^, The 
existence of glueballs is a natural prediction of QCD. The prediction for the glueball masses 
is, however, a non-trivial task. QCD sum rulesp, Q| and lattice QCdQ calculations obtain 
the lowest pseudoscalar glueball mass to be in the range of 1800 to 2600 MeV with lattice 
calculations giving a mass in the upper range. One cannot rule out the possibility that 
the resonant state X(1835) is a glueball based on our present understanding of the glueball 
masses alone. At this stage there is no compelling reasons to believe that the resonance 
is a pp bound state either |4j. These speculative particle states, although attractive, their 
existences are far from being established. More theoretical and experimental efforts are 
needed to go further. At a more modest level, even to know whether the data from BES 
can be consistently explained by a specific resonance and to further test the mechanism, 
more information about properties of the resonance is needed, such as how it is produced in 
radiative J/ ip decays and how it decays into other particles. 

In this work we study the properties of the X(1835) resonance using QCD anomaly and 
QCD sum rules assuming that this state is a pseudoscalar G p which couples strongly with 
two gluons which may or may not be a glueball or a pp bound state depending on whether 
this state has large mixing. We find that the matrix element < 0\GG\G p > is larger than 
< 0\GG\rj(rj') > indicating a large glue content in G p which is usually referred to as a glueball 



p 



in the literature. This leads to large branching ratios of (2.61 ~ 7.37) x 10 for J/ip — >• jG. 
and (2.21 ~ 10.61) x 10~ 2 for G p — > 7r + 7r — 77'. The coupling for G p — p — p interaction can 
also be determined. We finally discuss how other G p decay modes can be used to test the 
mechanism. 
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There have been considerable amount of literatures on production of a pseudoscalar in ra- 



diative J/i/j decays, in terms of QCD sum rulesj9| and perturbative QCD calculations |lll 111 | 



We follow the QCD sum rule approach in Ref . [9j such that we can treat radiative J / ip decays 
into rj, 1]' and G p in the same framework with QCD anomaly. In this framework, the radia- 
tive J/ip decay amplitudes are determined as follows: one first evaluates the internal charm 
quark loop contribution to the interaction of two-photon — > two-gluon, and then saturates 
the cc pair which couples to one of the external photons by J j ip and other resonant states 
using the standard procedure of the QCD sum rules. The two gluons are then converted into 
the related pseudoscalar states. This approach works the best when the final pseudoscalar 
has a mass-squared much smaller than Am 2 c ~ ^j/^,- For a pseudoscalar of mass 1833.7 MeV, 
there may be large corrections from significant two-photon to multi-gluon couplings since 
the factor ^G p / m j/V' ma J n °t be sufficient to suppress higher order contributions. However, 
one expects that the matrix elements of operators converting multi-gluon to a pseudoscalar 
G p must be smaller compared with that from the leading two gluon operator. The situation 
may not be too severe to damage the whole picture of the two gluon scenario. Also in our 
later discussions we will only use the ratios of two different J/ip — > jXi branching fractions, 
a large part of the uncertainty is expected to be cancelled out. One expects that the error 
range can be controlled to be within a factor of two. 

In this calculation the two-gluon operator with appropriate quantum numbers is, G^G^ . 
The matrix elements converting the two gluons into a pseudoscalar Xj is usually param- 
eterized as: fimf = (0|(3a s /47r)G / i iy G M!/ |X"i). Since the rest of the decay amplitude for 
J/ip — > 7X4 is independent of the final pseudoscalar state, the ratio of radiative branching 
fractions for X and Xj states is simply given by[l| 

BjJ/j, -> jJQ \f t m^ (1 ~ mym)^f 
13 2J(J/Vi-7*i) l/^ 2 | 2 (l-^ 2 M 2 W 3 ' U 

The parameters f v ,ri',G p play a crucial role in determination of J /if) — > •yGp in QCD sum 
rule approach. The parameters f VyV ' can be easily obtained from the QCD anomaly relations 
in the limit that the strange quark mass is much larger than the up and down quark masses. 
One hasQ 

(Ol^G^GHr?) = ^(cos^ - v^sinfl/oK, 

(OI^G^GHV) = y|( sin ^ + V2cos0/o)mJ, (2) 



where 9 is the rj-rj' mixing angle with rj = rjs cos 9 — rjo sin 9 and rj' = rjs sin 9 + 770 cos 9. fs,o 
are decay constants of the SU(3) octet rjs and singlet i] . 

There are many theoretical studies on the values of 9 and /g,o- 111 our study since only 
rj, 7]' and J/ip are involved, we will use related processes to determine /g o and 9. These 
processes include r\ — > 77, 7/ — > 77 and i?^/^ = B{J/ip — > ^yrj') / B{J /ip — > 777). We have 
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Using experimental values B(r] -> 77) = (39.43 ± 0.26)%, B(r/ -> 77) = (2.12 ± 0.14)%, 



B(J/ip -> 777) = (8.6 ± 0.8) x 10~ 4 and 5(J/V> 



777 



(4.31 ± 0.30) x 10- 3 |l3|, we 



obtain the ranges (central values) for the parameters as: 9 = —16.88° ~ — 18. 60°(— 17.72°), 
/ 8 = 0.98/, ~ 1.04/, (1.01/,) and / = 1.06/, ~ 1.21/, (1.08/,) with /, = 132 MeV being 
the pion decay constant. The correlations of these parameters are shown in Fig. 1. These 
values are consistent with the values determined from other considerations jl^ for 9 and /s,o- 
This gives us some confidence in using the QCD sum rule results for J/ip — > 777, J/ip — > 77/', 
and as well as for J/ip — > ^Gp. We will use the above values for 9 and /g o in our later 
discussions. 
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FIG. 1: The dependence of fo and /g on 9. The ranges are due to one a errors of data points. 

So far the parameter fa is not well understood. To obtain some information, we use the 
QCD sum rules to calculate it. The basic idea of QCD sum rule analysis in the present case is 
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to match the dispersion relation involving the hadronic spectral density p(s) to the vacuum 
topological susceptibility T(—q 2 ) —if d^xe iq ' x {Q\T[jj^x)ip S {Q)} |0) with the result found by 
using the operator product expansion. We follow Ref . |7J by making a Borel transformation on 
T(s) with J °° Im(T (s)e~ s / M2 ds / s and J °° Im(T(s))e~ s / M2 ds, and including the two ground 
pseudoscalar states r\ and rf, and G p in the resonant spectral density, to obtain the leading 
order matching conditions [7] 

fa y Gp e- mh * /M2 + fW^ mVM2 + f»- m *' IM2 

= I 1 bse- s / M2 ds + D 4 + 0(^-r)+inst., 
Jo M 

fc P < p e- m ^ /M2 + fK^ 12 + f>>- m2 ' IM2 

bs 2 e- s/M2 ds -D 6 + O(-r^) + inst., (4) 



where "inst" indicates direct instanton effects [1 51] . D± t % are related to the gluon condensa- 



tions, D 4 = 4(0|G^G^|0), and D 6 = 8g a f tt bc(0\Gl a G b l > a G c >'">\0), with D 4 = 7c{3a s /A7r) 2 D 4 , 
Dq = 7r(3a s /47r) 2 Z?6- b = (3a s /47r) 2 (2/7r)(l + 5a s /7c). In our numerical evaluations, we 
will use (0|« S G 2 |0) = (7.1 ± 0.9) x 1(T 2 GeV 4 and the relation (0\g 3 f abc G a G b G c \0) = 
1.2 GeV 2 (0|a s G 2 |0) Q- We comment that there are other 0"+ states around 1400 MeV 
region which may contribute to the spectrum density if these states contain large two gluon 
contents. We will assume that the gluon contents are small in these states and their contri- 
butions to the spectrum density can be neglected. 

To determine fa p , we take the usual practice to find the parameters fa p and s\ for a 
given Borel parameter M and look for a region where the dependence of fc p and s\ on M is 
insensitive. We will negelct the direct instanton effect in our calculation and will come back 
to comment on the effects later. Note that the analysis with a fixed G p mass here is different 
than previous ones where the mass of G p is taken as one of the parameters to be 

determined and therefore the results is in general different. The solutions for fc p depend on 
the value a s which we take to be the value at the scale /i = rriG p with a s = 0.35 ± 0.05. We 
find that solutions exist only for a restricted parameter space for a s and -D 4 6 . In certain 
ranges, for a given set of input values of a s and D^e, there are two solutions. For example 
with a s = 0.39, D 4 = 1.99 x 10~ 2 GeV 4 and As = 3.79 x 10~ 3 GeV 6 , we get the two solutions 
to be: a) s x = 3.2 GeV 2 and f Gp = 0.081 GeV, and b) si = 3.5 GeV 2 and f Gp = 0.091 GeV. 
When M is larger than 7 GeV, the solutions are fairly stable. As long as we choose an M 
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far above mc p , the power corrections from higher dimensional operators on the right hand 
sides of eq. (jH) can be suppressed. We note that for the solution with lower si, the value 
of Si is smaller than trq which cannot be considered to be a good solution since it implies 
that the continuum already starts to contribute to the sum rules in the resonant region 
in contradiction with the QCD sum rules assumptions. We therefore should choose the 
solution with the larger si. This solution allows a small gap between the resonances and the 
continuum. We show the results in Fig. 2 allowing a s to vary from 0.3626 (where solution 
begins to exist) to 0.4 (the one a allowed upper bound), and all other quantities, D^q, f v>ri i 
and 9 to vary within one a error ranges. We see that the dependence on M is very mild. 
We conclude that there are consistent solutions from QCD sum rules for a pseudoscalar of 
mass 1833.7 MeV, and obtain a conservative range for fa p with 

f Gp = 0.072 ~ 0.100 GeV. (5) 
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FIG. 2: fc p and si as functions of M. The upper and lower bounds are for one a ranges of D^q, 
fqtf and 6 with a s = 0.4 and a s = 0.3626, respectively. 
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When we obtained the range of fa„, the direct instanton effects were neglected. At 
large Q, the instanton effects are suppressed p, la bjj, for example spike distribution for 
instanton density results in an exponential suppression when Q 2 becomes larger than a GeV 2 
or so [6]. With the glueball mass fixed at 1835 MeV, the suppression may not be sufficient 
to neglect the contributions from the direct instanton. The effects of direct instanton may 
be substantial. However the detailed calculations depend on the instanton density and 
the density shape. A reliable evaluation of the instanton effects is difficult. Nevertheless 
model calculations show that instanton effects may be important due to modification to the 
normalization of the Wilson coefficients for the relevant operators. We will not go into the 
specific details as it is too model dependent. Some detailed discussions for direct instanton 
effects can be found in Ref.jy]. We emphasize that QCD sum rule results should be taken 
as an estimate within a factor of two. In our later discussions we will use the range for 
fa p obtained the above as a reference. Should a more precise value will be determined with 
some method, one can easily rescale the values accordingly. 

With the above range for fa we find that the matrix element < 0\a s GG\G p > is larger 
than < 0\a s GG\r](r]') > indicating that G p contains a large gluon content. With fa p deter- 
mined, we are now able to obtain information on the range for B(J/ip — > ^G p ) combining 
eq. (0) and experimental data on J/t|) -> 777(7/). We have 

B(J/ip -> 7 G P ) = (2.61 ~ 7.37) x 1CT 3 . (6) 

Using the BES data B(J/tp -> -fG p )B(G p -> tt+tt-t/) = (2.2 ± OA(stat) ± OA(syst)) x 
10 -4 , we can therefore also obtain the branching ratio of G p to 7r + 7r~r/ using our estimate 
of J/ijj — > •yGp. We have 

B(G P -> tt+tt-t]') = (2.21 ~ 10.61) x lCr 2 . (7) 

The branching ratio for this decay is large, but it does not contradict with known data. 

If the enhanced threshold pp events in J/?/; -> ^ypp is also due to the same G p state, we 
can obtain information about the interaction of G p with a proton and anti-proton pair, L = 
CgbP15PG p . Since the mass of G p is slightly below the threshold of two proton mass 2m p , 
one cannot simply take B(J/ip — > jG p •ypp) to be equal to B(J/ip — ► •yG p )B(G p — > pp). 
One must consider the off-shell effects of G p in terms of the Breit-Wigner approach. We 
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have 




(8) 



(1 - g 2 /m^)y y/j - ^y? 

Ami g (i-^H/^-my'+rG™^' 



In the above, we have assumed that Cgb does not depend on q 2 sensitively and has been 
moved out from the integration sign. 

Experimental data on B{J/ip — > jG p — > 'jpp) = (7.0 ± 0.4j| J) x 10~ 5 then implies 



One sees that threshold enhancement data in J/ip — > jpp can also be consistently explained. 

We now discuss how the resonance G p may decay into other particles. As have been 
pointed out earlier this state has a large gluon contents, it may be a glueball which is 
a SU(3) singlet. Of course one needs to be more open minded that it may has sizeable 
mixing with other state. We will take the state G p to be a glueball state and study the 
consequences from flavor symmetry point of view. One can easily study the branching 
ratios for J/ip — > ^Gp — > ^BB with BB a pair of octet baryons. The mixing effect can be 
easily implemented by introducing some mixing parameters. 

The coupling of a SU(3) singlet G p and octet baryon can be written as L = 
CgbG p TtB^B with SU(3) flavor symmetry. In Table 1 we list the ratios of riBB) = 
B(J/ip — > ^Gp — > -yBB)/B(J/ip — > jG p — > -ypp) for possible decay modes. We comment 
that the X(1835) contributions listed in Table 1 hold as long as the resonance is an SU(3) 
singlet and does not depend on the size of Cgb- If the resonance transforms non-trivially 
under the flavor SU(3) symmetry, the predictions would be different Q]. In principle experi- 
mental measurements of these branching ratios can provide important information about the 
nature of the resonance. However, the branching ratios for other baryon pair decay modes 
are much smaller than the branching ratio with a proton and anti-proton pair except the 
neutron and anti-neutron pair decay mode which is then experimentally difficult to carry 
out. The usefulness of these decay modes depends on whether, near the resonance region, 
the resonance contributions dominate over the non-resonance continuum parts which we will 
comment on later. 

A better test of the mechanism may come from other three meson decay modes of G p . 
To have further information on the branching ratios for G p decay into three meson modes, 



0.42 



0.85. 



(9) 
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1.58 x 1(T 3 


r(H-S-) 


1.38 x 1(T 3 


r(S+S+) 


1.45 x 1(T 2 


r(£°£°) 


1.38 x 1(T 2 


r(E~fr) 


1.28 x 1(T 2 


r(AA) 


4.23 x 1(T 2 


r(ren) 


0.96 



TABLE I: r(BB) for different J ftp -» 7GV -► decays. 

we here follow Ref. 19] to use U(3) chiral theory to describe how it couples to known meson 
particles. Notice that the use of U(3) symmetry will not change our previous discussions on 
G p — ► BB results. The reason we use U(3) chiral perturbation theory for the interaction is 
that it can naturally include many properties of chiral anomaly which our calculations for 
fi depend on. To the leading order there are four terms which may cause G p to decay [l^ 

L = iaxd^GpTr (E^E&E^E) 

+ ia 2 G p Tr ((E t 9 2 E — <9 2 E t E)<9 i ,E t <9 iy E) 
+ ia 3 G p Tr (xE - E f x) 

+ ia 4 G p Tr (x^d^d^ - xS^E^E) , (10) 

where E = exp[— i \/2M//] with / = f n /y/2 and M is the U(3) meson nonet, x is propor- 
tional to the light quark masses and is given by x — diag{m Q K1 m%, 2m 2 K — m 2 ). The last two 
terms in eq. (jTUjl comes from explicit U(3) (and SU(3)) breaking due to quark masses. 

The two SU(3) breaking terms in eq. (jl(Jj) . if dominant, will lead to the main decay mode 
to be G p -»■ TrfsTi^Q with a very suppressed rate for G p — > 7i + n r}' . The BES data indicates 
that the decay mode G p — > 7i + 7i^r] f has a large branching ratio compared with other three 
particle decays (yet to be discovered), therefore, these two terms may be suppressed. If 
the first two terms dominate, we obtain the effective Lagrangian for the decay amplitude 
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G p — > 7r + 7r 7]' to be, 

L = -^(V2cos9 + sm9) (a 1 d l ,G p (d'"r]'d l/ n + d 1/ TT- 
V3/ 3 

+ d»it-d v Tt + d v ri' + d"TT + d v r ] 'd I/ TT-) 

+ a 2 G p {d 2 r(d v ix+d u ix~ + d 2 it-d v Tt + d v r]' 

+ df 2 n + d v ri'd u n-)) . (11) 
which leads to a decay amplitude, 

+ _ , y^cosfl + sinfl 

M (G p - 7T+7T 77) = 

x (ai[(s - m 2 ^ - m 2 v ,)(s - 2m 2 ) 
+ (* - m G P - m l)(t ~ m l' ~ m l) 

+ ( u - m G p - m l)( u - m J' - m l)] 

+2a 2 [m 2 ,(s - 2ml) + m ^(* ~ m J/ - m ^) 

+m 2 r ( M -m 2 / -771 2 r )]), (12) 

where s = + P7,--) 2 , * = (p n - + Pr, 1 ) 2 an d M = (Pn+ + Pr, 1 ) 2 ■ 

To obtain the BES result for J/ip — > 77r + 7r~7/ with the B(J/ip — > 7G P ) obtained earlier, 
we have 

|«i | 2 - 0.376i?e(aia;) + 0.038|a 2 | 2 = (0.94 ~ 4.52) x 10~ 5 GeV~ 2 . (13) 

Since there are two parameters, ai and a 2 , involved, with data from J/xjj — > 7G P — > 
77r + 7r~7/ alone it is not possible to predict branching ratios for other three meson decay 
modes. There are several kinematically allowed decay modes which can provide more 
information. These additional decay modes include: r]'7T°n , r]7i + 7Y^, r]n 7r°, r/K + K~, 
r]K°K°, ir°K°K°, ir°K+K-, tt + K°K-, ir~K°K+, and 777777. With the assumption that 
the dominant contributions come from the 01,2 terms, even with SU(3) breaking ef- 
fects from the final meson mass differences, one would have the following relations be- 
tween branching ratios: B(7r 7r°r)'(i])) = B(7i + n~r)'(r)))/2, B(K + K~i]) = B(K°K i]), and 
B(n°K + K~) = B(ir°K°K ) = B(n'K + K°)/2 = B(n + K~K°)/2. We list the ratios 
r{piP2P3) = B(G P — > P1P2P3) I 'B{G P — > ir + ir~r)') in Table 2. The branching ratio for 
G p — > 7r + 7r _ 77 can be larger than that for G p — > ir + ir~r)'. No observation of J/ip — > 77r + 7r~77 
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near the resonance would imply that there may be large cancellations between terms pro- 
portional to di and a 2 for this decay mode. G p — > 7rKK are some another modes which 
may have large branching ratios near the resonance. The branching ratios here are different 
from those predicted if the resonance transforms non-trivially under the SU(3) symmetry^]. 
These new decay modes can serve to test the mechanism proposed in this paper. We urge 
the BES collaboration to carry out an systematic analysis to obtain more information about 
the properties of the resonant state X(1835). 



r(ir + ir rj) 


11.26(|ai| 2 - 0.186i?e(ai4) + 0.009|a 2 | 2 )/6(7r+7r~y) 


r(K + K-rj) 


1.57 x 10~ 3 (|ai| 2 - 0A74Re(a ia * 2 ) + 0.056|a 2 | 2 )/6(7r+7r-?/) 


r(K + K~-K°) 


2.47(|ai| 2 -0.291^e(aia^) +0.021|a 2 | 2 )/6(vr + 7r-?/) 


r (777777) 


0.006(|ai| 2 - 0.538iie(aia^) + 0.073|a 2 | 2 )/6(7r + 7r-?/) 



TABLE II: r{pip2p-z) for various G p — > P1P2P3 decay modes. Here b(ir + iT rf) = \a± 
0.376Re{cna* 2 ) + 0.038|a 2 | 2 . 



In the above discussions about predictions for other decay modes, only contributions from 
the resonance are included. There are also non-resonance effects. To isolate the resonance 
contributions near the resonance region, one should remove the off-resonance contributions 
by extrapolating the data away the resonance to the resonance region. Theoretical study 
of the non-resonance contributions is more complicated. We briefly discuss how the non- 
resonance contributions for J/ip —> ^BB and J/ip — > ■77r + 7r — 77(77') can be parameterized 
using SU(3) flavor symmetry. The general form for non-resonance J/ip — > 7I? B has been 
discussed in Ref . [ls| . As far as the SU(3) structure is concerned, that is, neglecting Lorentz 
structure, we have the amplitudes for non-resonance contributions to J/ip — > 'jBB and 
J/ij) — > P1P2P3 to be given by 

M(J/?P -> 1 BB) ~ Tr[B(D 7 {Q, B} + F 7 [Q, B})), 

M(J/tP -> 7P1P2P3) ~ arnTr^QM 2 ) + +bi] 1 r] 1 Tr(QM) + cTr(QM 3 ). (14) 

where Q = diag(2/3, —1/3, —1/3) is the electric charge matrix. D 1 , F 7 , a, b and c are form 
factors. 

The complications in extracting information on the form factors _D 7 and F 7 , and, a, b and 
c are two folds. One of them is that the Lorentz structure is much more complicated, that is 
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for each of the form factors _D 7 , _F 7 , a and 6, there are actually several of them depending on 
how the Dirac matrices are inserted in the BB bi-baryon product, and how derivatives are 
taken on the meson and baryon fields. Another is that the form factors in general depend 
on combinations of the invariant masses of pairs of particles in the final state. Nevertheless, 
fitting data, one can obtain information about the parameters. Measurements of branching 
ratios of J/ip — > / -fBB, / ypip2P3 via non-resonance can also provide important information 
for understanding the whole physics picture. 

In summary, we have studied some implications of a possible h resonance in the ir + 7i^r]' 
spectrum in J/ip — > 77r + 7r~7/ observed by BES. We have shown that it has a sizeable matrix 
element for < 0\GG\G p >. The branching ratios for J/ip — > '-/Gp and G p — > 7r + 7r~?/, 
using QCD anomaly and QCD sum rules, are determined to be (2.61 ~ 7.37) x 10~ 3 and 
(2.21 ~ 10.61) x 10~ 2 , respectively. The coupling for G p —p—p interaction is also determined. 
We conclude that a pseudoscalar 0~ + with large gluon content can consistently explain 
the data. We have also studied branching ratios of other decay modes using SU(3) flavor 
symmetry. We find that J/xjj — ► jG p — > 7(7r + 7r~?7, KKn°) can provide useful tests for the 
mechanism proposed here. 
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